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Abstract 

The AMS-02 collaboration has released preliminary data on the antiproton fraction 
in cosmic rays. The surprisingly hard antiproton spectrum at high rigidity has triggered 
speculations about a possible primary antiproton component originating from dark matter 
annihilations. In this note, we employ newly available AMS-02 boron to carbon data to 
update the secondary antiproton flux within the standard two-zone diffusion model. The 
new background permits a considerably better fit to the measured antiproton fraction 
compared to previous estimates. This is mainly a consequence of the smaller slope of the 
diffusion coefficient favored by the new AMS-02 boron to carbon data. 



1 Introduction 


The AMS-02 experiment at the International Space Station is performing a high precision 
measurement of charged cosmic ray fluxes. The spectrum of antiprotons is of particular interest 
as it provides a very sensitive probe of possible exotic sources in the galaxy. In this light, 
preliminary pj-p data [1] which indicate a nearly constant ratio at high rigidity have caused 
considerable excitement. Interpretations of the pjp data in terms of TeV mass dark matter 
annihilations have been suggested [2-5].^ 

On the other hand, the p flux suffers from sizeable uncertainties related e.g. to the modeling 
of cosmic ray transport in the galaxy. It was pointed out in [9,10] that after fully accounting for 
uncertainties, no unambiguous excess is left in the AMS-02 pjp ratio. Still, viable configurations 
consistent with the data appeared somewhat at the edge of the uncertainty band provided in [9] . 

In this note we improve the calculation of the secondary p background^ within the two- 
zone diffusion model. This is achieved by employing the new AMS-02 boron to carbon (B/C) 
data [11]. B/C being a secondary to primary ratio [12] can be used to estimate propagation 
parameters of the intergalactic transport equation. The AMS-02 data considerably extend the 
rigidity range compared to previous measurements and allow for a more robust determination 
of propagation parameters which - at high energies - are less affected by degeneracies. Our 
approach is more data-driven compared to previous work as we do not attempt to predict the 
fluxes of primary cosmic ray nuclei within the diffusion model. Merely, we directly take the 
observed primary fluxes as an input for the calculation of secondary cosmic ray production. 
We also make use of the updated p production cross sections [13] which were recently derived 
from scattering data of the NA49 experiment [14,15] . 

The secondary antiproton flux obtained with the new sets of propagation parameters tends 
to decrease less strongly with energy compared to previous estimates [9,13,16,17]. This is 
mainly triggered by the smaller power law index of spatial diffusion b favored by the AMS-02 
B/C data. Intriguingly, the allowed range of b includes the theoretically favored values b = 0.33 
and 5 = 0.5 which correspond to a Kolmogorov and Kraichnan type power spectrum of the 
galactic magnetic field. 

The harder antiproton spectrum following from our propagation analysis considerably im¬ 
proves the consistency of the AMS-02 p/p data with the background hypothesis. 

2 Transport of Cosmic Rays 

The propagation of all cosmic ray species can be performed in a unified framework. It is 
common to start with a parameterization of primary sources and then to calculate the fluxes 
of all primary and secondary species resulting from the network of scattering and spallation 
reactions in the galactic disc [18-21]. As we are mainly concerned with the secondary cosmic 
ray species B and p, we instead perform a more data-driven approach and use observed primary 
fluxes as an input for our calculation. Before discussing this in more detail, we will very briefly 
review the relevant aspects of cosmic ray propagation. 

^See also [6-8] for recent dark matter studies using antiprotons. 

^Secondary cosmic rays are generated by the scattering of primary cosmic rays on the interstellar matter. 
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2.1 Propagation in the Diffusion Model 

Cosmic rays propagating through the galaxy are affected by various processes which are en¬ 
coded in the diffusion equation. The latter relates the production rate (source term) qi of the 
cosmic ray species i to the resulting space-energy density iVj 

V(-K VlVi + Vc Ni) + dEihot Ni - KeE dENi) + r^nn N, = q,, (1) 

where E stands for the total energy. K accounts for diffusion on magnetic field inhomogeneities. 
The galactic wind Vc is responsible for convection. The function 6 tot includes Coulomb, 
ionization and adiabatic energy losses as well as energy losses induced by reacceleration. Kee ^^ 
the energy diffusion coefficient and Tann the annihilation rate on the interstellar matter [22,23]. 

In this work we employ the two-zone diffusion model [16, 18] in which diffusion occurs 
homogeneously in a cylinder of half-height L around the galactic disc.^ K is parameterized as 

where Kq and 5 are the normalization and the power law index of spatial diffusion. Magne- 
tohydrodynamics considerations suggest values 5 = 0.33 (Kolmogorov spectrum) or 5 = 0.5 
(Kraichnan spectrum) depending on the detailed modeling of magnetic field turbulences (see 
e.g. [24]). The velocity and rigidity of the cosmic ray particle are denoted by /3 and TZ. The 
convective wind Vc is taken to be constant in the diffusion cylinder and pointing away from the 
galactic disc. Finally, the Alfven speed of magnetic shock waves 14 enters Kee and completes 
the set of five independent propagation parameters: L, Kq, 5, 14 and 14- 


2.2 Source Terms 


A common approach to predict cosmic ray spectra is to start from a parameterization of 
primary sources. A standard choice is oc with t/ = 0 [18] or = 1 [21] (see 

also [25]). At high energies, this leads to a power law form of the primary flux oc . The 
index 7 is matched to observed high energy fluxes, e.g. 7 = 2.65 — <5 in [21]. Primary fluxes 
over the full energy range are then calculated by means of the diffusion equation and seeded 
into the source terms for secondaries. 

One drawback of this approach is that the low energy parameterization of primary sources 
is rather ambiguous. The observed spectral breaks in the proton and helium spectra [26, 27] 
further indicate that a simple power law may not be sufficient to model the primary sources 
at high energies. For this reason, we have chosen not to start from a parameterization of 
primary sources, but directly use the observed primary fluxes as an input for the calculation 
of secondaries. 

Secondary cosmic rays like B or p are generated by scattering of (mainly) primary cosmic 
rays on the interstellar matter. The source term for secondary production reads 


vTiT) = Y. 

j 


dT' 


nu 



+ nue 





(3) 


Here <h^®(r') is the interstellar flux of the incoming cosmic ray species j as a function of its 
kinetic energy per nucleon T', while T denotes the kinetic energy per nucleon of the outgoing 

^See [19,20,22] for alternative approaches to solve the diffusion equation. 
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secondary particle i. The differential cross section da/dT for the production of i includes 
the prompt process as well as the production via intermediate unstable particles which decay 
into i. We consider only scattering processes off hydrogen and helium in the galactic disc 
and use constant values nn = 0.9 cm“^ and nue = 0.1 cm“^ for their number densities [21]. 
Furthermore, we take to be radially constant in the disc. Local secondary cosmic ray fluxes 
are not particularly sensitive to the last two assumptions as long as they are used in the same 
way for the determination of propagation parameters. 


2.3 Fluxes of Primary Cosmic Rays 

The interstellar fluxes of primary cosmic rays have been measured over a wide energy range 
by various experiments [28]. As experimental data refer to the top-of-the-atmosphere (TOA) 
fluxes we have to take into account solar modulation. We ignore charge-sign dependent effects 
and use the force field approximation [29]. The value of the force field 4> for experimental 
data prior to 2011 is extracted from [30]. The extended table [31] further allows us to deduce 
(f) = 0.57 GV for the new AMS-02 data. 

The proton and helium fluxes released by the AMS-02 collaboration [26, 27] clearly in¬ 
dicate a spectral hardening at rigidity TZ ~ 400 GV consistent with previous measurements 
by PAMELA [32] and CREAM [33] (see also [34-37]). We thus parameterize the interstellar 
proton and helium fluxes as the sum of two power laws 

^P(r)=(^) (AiT-^i+AsT-^^). (4) 

The bracket which multiplies the two terms accounts for the low energy (T < 10 GeV/n) 
behavior of the observed fluxes. We determine the free parameters by a fit to the data of 
AMS-02 and CREAM. The best fit values are presented in table 1. The experimental data are 
shown with our fit function in figure 1 . 

Ai A 2 b c 71 72 

Protons 53488 3294 2A0 2.588 3.168 2.576 

Helium 3194 141.5 5.18 1.482 3.163 2.451 

Table 1: Fit parameters for the primary spectra of protons and helium (4) with Ai,A 2 in 
m“^s“^sr“^(GeV/n)~^ and b in GeV/n. 


For the secondary production of boron the most relevant process is the spallation of carbon 
in the interstellar medium. However, several other elements, namely oxygen, nitrogen, neon, 
magnesium and silicon contribute significantly [23] . The primary fluxes are determined from a 
fit to the experimental data below lOOOGeV/n. In order to cover the low, intermediate and high 
energy range, we extract the measurements of ACE [38], HE AO [39] and CREAM-H [40] from 
the cosmic ray database [28]. Eor the carbon flux an additional data set from PAMELA [41] is 
considered. Primaries above lOOOGeV/n are excluded from our fit as they do not contribute to 
the boron flux in the energy range accessible to AMS-02. The interstellar flux as a function of 
the kinetic energy per nucleon T is modeled as in the case of protons and helium, but without 
the second power law (as the presence of a spectral break is not yet statistically significant) 


4'P(r) = A 


T + b 


C 

T-'y 


(5) 
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Figure 1: Proton and helium fluxes measured by AMS-02 [26,27] and CREAM [33]. Also 
shown are the fit functions employed in this work. The latter were modulated with a force 
held (/) = 0.57 GV. 


In order to avoid a too strong impact of low energy data points, suffering from uncertainties 
like solar modulation, we first determine the power law index from the data above 10 GeV/n. 
In the second step, we fit the remaining parameters to the full data set. The best fit parameters 
for oxygen, nitrogen, carbon, neon, magnesium and silicon are presented in table 2 and the 
corresponding fluxes are displayed in figure 2. 



A 

b 

c 

7 

Carbon 

23.13 

1.12 

2.086 

2.666 

Oxygen 

30.90 

1.76 

1.844 

2.748 

Nitrogen 

15.67 

1.39 

2.759 

3.077 

Neon 

5.01 

1.57 

2.035 

2.769 

Magnesium 

7.72 

1.97 

1.927 

2.829 

Silicon 

5.33 

2.04 

1.822 

2.746 


Table 2; Fit parameters for the primary spectra described by (5) with A in 
m“^s“^sr“^(GeV/n)~^ and b in GeV/n. 


3 Propagation Parameters from a B/C Analysis 

In the next step, we wish to determine sets of propagation parameters consistent with the 
preliminary B/C spectrum released by AMS-02 [11]. For the calculation of the and ^®B 
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Figure 2: Measured primary cosmic ray fluxes. Also shown are the fit functions employed in 
this work. The latter were modulated with the force held 4> = 0.43 GV for ACE (solid line) 
and the force held (f) = 0.74 GV for HEAO (dashed line). 
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source terms we consider spallation of the cosmic ray species listed in table 2 which account 
for ~ 98% of the boron production [23]. The isotopic composition of the incoming cosmic rays 
is extracted from the cosmic ray data base [28] . 

For the nuclear spallation cross sections, we use the straight-ahead approximation which 
implies that the kinetic energy per nucleon is preserved during a spallation process. Cross 
sections of nuclei on hydrogen are taken from [42] (as provided in the DRAGON code [20]), 
the energy-dependent enhancement factor for spallation on helium is extracted from [43]. We 
include the indirect production of boron via intermediate radioactive isotopes in the spallation 
cross sections and also account for the production of via neutron stripping of For the 
annihilation cross sections on the interstellar hydrogen and helium we use [44,45]. 

Spallation cross sections suffer from uncertainties related to the lack of experimental data 
at high collision energies (see [46] for a comprehensive discussion). In order to model this 
uncertainty, we follow [47] and allow for a systematic energy bias in the parameterization of 
spallation cross sections. This is done by multiplying cross sections by a factor [r/(GeV/n)]^ 
at T > 1 GeV/n which goes into a constant at T = 10 GeV/n. As can be seen in figure 3 of [47] 
any choice between A = —0.05 and A = 0.05 is consistent with existing experimental data. 

For our analysis, we generate a random set of propagation parameters within the intervals 
L = 2 — 15 kpc, 6 = 0.2 — 0.9, 14 = 0 — 20 km s“^, 14 = 0 — 100 km s“^ and free Kq. For the 
given conhguration we determine the B/C ratio after accounting for solar modulation. Then, 
we perform a t^st against the B/C data measured by AMS-02 [11]. The total experimental 
uncertainty is extracted from slide 14 of [11] by adding errors in quadrature. In order to 
account for the uncertainty in the spallation cross section, we allow the conhguration to choose 
any value between A = —0.05 and A = 0.05. A conhguration is selected if y^/d.o.f. < 2 for 
the optimal A. 

In hgure 3 we depict the envelope of B/C ratios within a sample of 500 selected con- 
hgurations. The set of propagation parameters which yields the best ht among this sample 
(x^/d.o.f. = 1.2) is given in table 3. 

6 Kq (kpc^ Myr“^) L (kpc) 14 (km s“^) 14 (km s“^) 

0.408 0.0967 ITT OA 31)9 

Table 3: Propagation parameters giving the best ht to the AMS-02 B/C data among a sample 
of 500 selected conhgurations (see text). 

As the presence of convective winds in the galaxy is speculative [24] it is remarkable that the 
best conhguration has almost vanishing convective wind. This already implies that the AMS-02 
B/C data are consistent with a pure diffusion-reacceleration model of cosmic ray propagation. 
The value of 6 is of particular interest as it sets the slope of secondary to primary ratios at high 
rigidity. In hgure 5 we show 6 and 14 for our sample of conhgurations. There is an obvious 
correlation between the two parameters: larger values of <5 require a higher convective wind. 
Despite the scatter and the difference in the statistical method, we can make out a clear trend 
if we compare our results with previous B/C analyses in the two-zone diffusion model [18,21]. 
The allowed parameter range lies at smaller 5. In our sample we hnd a median of <5 = 0.48, 
while the best ht conhguration of [21] has 6 = 0.86. The mentioned analysis is, however, 
based on earlier B/C data sets. We have checked that the best ht conhgurations of [18,21] 
substantially underestimate the AMS-02 B/C ratio at high rigidity. As cosmic ray transport at 
high energies is dominated by diffusion, the B/C ratio in this regime is particularly important 
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Figure 3: Envelope of B/C ratios within the sample of 500 selected configurations (see text) 
together with the experimental data of AMS-02. The B/C ratio for the best fit configuration 
in the sample (table 3) is indicated by the dotted line. 


for a robust determination of 6. The good accuracy of the AMS-02 data at high rigidity gives 
us confidence that the shift towards lower 6 will be confirmed by future data sets. This is 
further substantiated as analyses of the PAMELA [41] and earlier AMS-02 B/C data sets [46] 
which were performed in slightly different propagation models yield a range of S compatible 
with our findings. 

The range of 6 found in our analysis includes the theoretically favored values 6 = 0.33 
and S = 0.5. We will show in the next section that the lower S does also have important 
implications for the antiproton fraction. 

4 Antiproton Fraction 

Secondary antiprotons originate from the scattering of protons and helium on the interstellar 
matter. The primary fluxes are given in (4) and the fit parameters in table 1. Contributions 
from heavier nuclei are negligible. 

We make use of the new calculation of antiproton production cross sections performed in [13] 
(for other recent approaches see [48,49]). Compared to previous parameterizations [50,51] it 
contains a detailed treatment of antiproton production via hyperon decay as well as possible 
isospin effects in antineutron production. Eurther the description of processes involving helium 
has been improved. We follow [52] and include tertiary antiproton production by inelastic 
scattering of secondary antiprotons on the interstellar matter with annihilation and inelastic 
cross sections taken from [53,54]. 

Eor the sample of 500 configurations consistent with B/C (see previous section) we de¬ 
termine the interstellar p flux. We take the envelope of all resulting fluxes to model the 
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Figure 4: Antiproton fraction predicted from pure secondary production compared to the 
AMS-02 data. The inner band encompasses propagation uncertainties (see text), the full band 
also includes uncertainties in the p production cross sections. The antiproton fraction for the 
propagation configuration (within our sample) which yields the best fit to B/C (table 3) and 
for the configuration which yields the best fit to the p/p data are indicated by the dotted and 
the dashed line, respectively. 


propagation uncertainties. The corresponding p/p ratio after accounting for solar modulation 
{(j) = 0.57 GV) is shown with the AMS-02 data in figure 4. The broader band in the same 
figure is obtained by including the uncertainties in the antiproton production cross sections 
from [13]. The p/p ratio for the configuration of table 3 is also shown. 

It can be seen that the secondary antiproton background is in good agreement with the 
data, primary sources of antiprotons are not favored. To make this more explicit, we have 
performed a test against the AMS-02 pjp data for each configuration within our sample. 
Even before taking into account the uncertainties in the antiproton production cross section, 
we find a configuration with y^/d.o.f. as low as 0.5. The p/p ratio for this configuration is also 
shown in figure 4. 

In figure 5, which shows 5 and V/ for the sample of configurations selected in the B/C 
analysis, we have marked those which are also consistent with the p/p data. As a criterion we 
again required^ x^/d.o.f. < 2. 

There is a trend that the AMS-02 p/p data favor those sets of propagation parameters with 
smaller 6. This is a consequence of the almost flat shape of the AMS-02 p/p ratio at high 

^We also took into account the uncertainty in the production cross section. For each set of propagation 
parameters, we calculated the minimal, medium and maximal flux within the cross section uncertainty band. 
Then, we defined a parameter which smoothly interpolates between the three fluxes and selected the parameter 
which minimizes If for this optimal choice of the production cross section y^/d.o.f. < 2 the configuration is 
taken to be consistent with the p/p data of AMS-02. 
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Figure 5: Colored points indicate values of 6 and Vc within a sample of 500 configurations 
consistent with the AMS-02 B/C data. Configurations which simultaneously fit the AMS-02 
p/p data are displayed in red. Theoretically preferred values of 5 = 0.33 and 6 = 0.5 are 
indicated by the dashed lines. 


consistent with B/C 
consistent with B/C and p/p 



rigidity. Both AMS-02 data sets, B/C and p/p, are consistent with Kolmogorov (5 = 0.33) 
and with Kraichnan (<5 = 0.5) type diffusion. In the case of d = 0.33 convective winds are not 
required. 

While we have shown that the AMS-02 p/p data are consistent with the background hy¬ 
pothesis, the relatively large uncertainties still leave considerable room for a possible primary 
p component. In order to improve the search for dark matter in cosmic rays, it will be crucial 
to shrink the error bands in the background. Among other things, this requires new experi¬ 
mental data on nuclear spallation cross sections which are a limiting factor in determining the 
propagation parameters. 

5 Conclusion 

We investigated whether the preliminary AMS-02 p/p data can be explained by the secondary 
antiproton background within the two-zone diffusion model of cosmic ray propagation. In the 
first step we determined viable sets of propagation parameters by use of the new AMS-02 B/C 
data. Compared to previous analyses, we found a trend towards smaller slopes of the diffusion 
index. The theoretically favored values 5 = 0.33 and 5 = 0.5 are both consistent with the data. 

We then determined the error band of the p/p ratio by including uncertainties in the 
propagation parameters as well as in the cross sections for antiproton production. The AMS- 
02 p/p data reside well within the error band. No primary sources of antiprotons like dark 
matter annihilations are required. Indeed, the background we obtained gives a somewhat better 
fit to the AMS-02 data compared to the background of [9]. This is mainly a consequence of 
the smaller 5 favored by our propagation analysis, which translates into a harder antiproton 
spectrum at high rigidity. 
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AMS-02 is expected to further extend the energy range and the accuracy of the antiproton 
data within the next years. The secondary antiproton ratio and the uncertainty band provided 
in this work will allow an improved background modeling in future dark matter studies. 

Upcoming AMS-02 data on primary and secondary fluxes will further reduce the uncertain¬ 
ties in cosmic ray propagation. In the near future, uncertainties in the cross sections entering 
secondary source terms will be a limiting factor. This includes nuclear spallation cross sec¬ 
tions which are not well-measured at high energies as well as the antiproton production cross 
sections (in particular the component from antineutron decay). Future dark matter searches 
in cosmic rays will aim at identifying a sub-dominant signal on top of a large astrophysical 
background. In this light, there is an urgent need for a dedicated experimental program [55] 
to measure cross sections in the relevant energy range. 
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